Quantum Spin Hall Effect in Graphene 
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We study the effects of spin orbit interactions on the low energy electronic structure of a sin- 
gle plane of graphene. We find that in an experimentally accessible low temperature regime the 
symmetry allowed spin orbit potential converts graphene from an ideal two dimensional semimetal- 
lic state to a quantum spin Hall insulator. This novel electronic state of matter is gapped in the 
bulk and supports the transport of spin and charge in gapless edge states that propagate at the 
sample boundaries. The edge states are non chiral, but they are insensitive to disorder because 
their directionality is correlated with spin. The spin and charge conductances in these edge states 
are calculated and the effects of temperature, chemical potential, Rashba coupling, disorder and 
symmetry breaking fields are discussed. 

PACS numbers: 73.43.-f, 72.25.Hg, 73.61. Wp, 85.75.-d 
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The generation of spin currents solid state systems has 
been a focus of intense recent interest. It has been argued 
that in doped semiconductors the spin orbit (SO) inter- 
action leads to a spin- Hall effect U 0, in which a spin 
current flows perpendicular to an applied electric field. 
The spin Hall effect has been observed in GaAsHli. Mu- 
rakami et al. have identified a class of cubic materials 
that are insulators, but nonetheless exhibit a finite spin 
Hall conductivity. Such a "spin Hall insulator" would be 
of intrinsic interest, since it would allow for spin currents 
to be generated without dissipation. 

In this paper we show that at sufficiently low energy 
a single plane of graphene exhibits a quantum spin-Hall 
(QSH) effect with an energy gap that is generated by 
the SO interaction. Our motivation is twofold. First, 
Novoselov et al. @ have recently reported progress in the 
preparation of single layer graphene films. These films 
exhibit the expected ambipolar behavior when gated and 
have mobilities up to 10 4 cm 2 /Vs. Thus, the detailed ex- 
perimental study of graphene now appears feasible. We 
believe the QSH effect in graphene is observable below a 
low but experimentally accessible temperature. Secondly, 
we will show the QSH effect in graphene is different from 
the spin hall effects studied for three dimensional cubic 
systems in Ref. because it leads to a phase which is 
topologically distinct from a band insulator. The QSH 
effect in graphene resembles the charge quantum Hall ef- 
fect, and we will show that spin and charge currents can 
be transported in gapless edge states. As a model sys- 
tem, graphene thus identifies a new class of spin Hall 
insulator. It may provide a starting point for the search 
for other spin-Hall insulators in two dimensional or in 
layered materials with stronger SO interaction. 

SO effects in graphite have been known for over 40 
years 0, and play a role in the formation of minority 
hole pockets in the graphite Fermi surface 8J. However, 
these effects have largely been ignored because they are 
predicted to be quite small and they are overwhelmed by 
the larger effect of coupling between the graphene planes. 
Unlike graphite which has a finite Fermi surface, however, 



graphene is in a critical electronic state which can be 
strongly affected by small perturbations at low energy. 

Graphene consists of a honeycomb lattice of carbon 
atoms with two sublattices. The states near the Fermi 
energy are 7r orbitals residing near the K and K 1 points 
at opposite corners of the hexagonal Brillouin zone. An 
effective mass model can be developed 9] by writing the 
low energy electronic wavefunctions as 

*(r) = [{u ak ,u B k), (uak>,ubk>)} ■ ip(r) (1) 

where u^ B)(K.K')( r ) describe basis states at momentum 
k = K, K' centered on atoms of the A, B sublattice. tp(r) 
is a four component slowly varying envelope function. 
The effective mass Hamiltonian then takes the form, 



Ho 



-ihv F ip\cr x T z d x + o-ydy)ip. 



(2) 



Here <x and f are Pauli matrices with a z — ±1 describing 
states on the A(B) sublattice and t z = ±1 describing 
states at the K(K') points. This Hamiltonian describes 
gapless states with E(q) — ±i!_p|q|. Without spin, the 
degeneracy at q = is protected by symmetry. The only 
possible terms that could be added to open a gap are pro- 
portional to a z or a z r z . The a z term, which corresponds 
to a staggered sublattice potential is odd under parity 
(which interchanges the A and B sublattices). The (t z t z 
term is even under parity, but odd under time reversal 
(which interchanges K and K'). 

The SO interaction allows for a new term, which will 
be the focus of this paper: 



Hso = A so ^a z T z s z ip. 



(3) 



Here s z is a Pauli matrix representing the electron's spin. 
This term respects all of the symmetries of graphene, and 
will be present. Below we will explicitly construct this 
term from the microscopic SO interaction and estimate 
its magnitude. If the mirror symmetry about the plane 
is preserved then this is the only allowed spin dependent 
term at q = 0. If the mirror symmetry is broken (either 
by a perpendicular electric field or by interaction with a 



2 



substrate) then a Rashba term[lO| of the form (s x p) ■ z 
is allowed, 

'Hr = XR^(o- x r z Sy - cr y s x )i>. (4) 

For \ji = 0, A so leads to an energy gap 2A so with 
E(q) = ±y / (hv F q) 2 + A 2 so . For < X R < A so the en- 
ergy gap 2(A so — Xr) remains finite. For Xr > A so the 
gap closes, and the electronic structure is that of a zero 
gap semiconductor with quadradically dispersing bands. 
In the following we will assume that Xr < A so and an- 
alyze the properties of the resulting gapped phase. This 
assumption is justified by numerical estimates given at 
the end of the paper. 

The gap generated by o~ z t z s z is different from the gap 
that would be generated by the staggered sublattice po- 
tentials, <j z or a z s z . The ground states in the presence 
of the latter terms are adiabatically connected to sim- 
ple insulating phases at strong coupling where the two 
sublattices are decoupled. In contrast, the gap param- 
eter a z T z s z produces gaps with opposite signs at the K 
and K' points. This has no simple strong coupling limit. 
To connect smoothly between the states generated by a z 
and o~ z t z s z one must pass through a critical point where 
the gap vanishes, separating ground states with distinct 
topological orders. 

The interaction (3) is related to a model introduced 
by Haldane[iH as a realization of the parity anomaly in 
(2+1) dimensional relativistic field theory. Taken sepa- 
rately, the Hamiltonians for the s z = ±1 spins violate 
time reversal symmetry and are equivalent to Haldane's 
model for spinless electrons, which could be realized by 
introducing a periodic magnetic field with no net flux. 
As Haldane showed, this gives rise to a g z t z gap, which 
has opposite signs at the K and K 1 points. At tempera- 
tures well below the energy gap this leads to a quantized 
Hall conductance a xy — ±e 2 /h. This Hall conductance 
computed by the Kubo formula can be interpreted as the 
topological Chcrn number induced by the Berry's curva- 
ture in momentum space [T3. Ilijj . Since the signs of the 
gaps in (3) are opposite for opposite spins, an electric 
field will induce opposite currents for the opposite spins, 
leading to a spin current 3 S = (?i/2e)(J| — JjJ charac- 
terized by a quantized spin Hall conductivity 



Since spin currents do not couple to experimental probes 
it is difficult to directly measure (5). Moreover, the con- 
servation of s z will be violated by the Rashba term (4) as 
well as terms which couple the ir and a orbitals. Nonethe- 
less, Murakami et al. [14| have defined a conserved spin 
s 2 ( c ), allowing <j s xy to be computed via the Kubo formula. 
We find that a% y computed in this way is not quantized 
when Xr ^ 0, though the correction to (5) is small due 
to carbon's weak SO interaction. 



In the quantum Hall effect the bulk topological order 
requires the presence of gapless edge states. We now show 
that gapless edge states are also present in graphene. We 
will begin by establishing the edge states for Xr = 0. 
We will then argue that the gapless edge states persist 
even when Xr ^ 0, and that they are robust against 
weak electron electron interactions and disorder. Thus, 
in spite of the violation of (5) the gapless edge states 
characterize a state which is distinct from an ordinary 
insulator. This QSH state is different from the insulators 
discussed in Ref. yj which do not have edge states. It 
is also distinct from the spin Hall effect in doped GaAs, 
which does not have an energy gap. 

For Xr = 0, the Hamiltonian (2,3) conserves s z , 
and the gap less edge states follow from Laughlin's 
argument |l5|. Consider a large cylinder (larger than 
Kvf/ 'A so ) and adiabatically insert a quantum (f> — h/e of 
magnetic flux quantum down the cylinder (slower than 
A so /?i). The resulting azimuthal Faraday electric field 
induces a spin current such that spin fi is transported 
from one end of the cylinder to the other. Since an adia- 
batic change in the magnetic field cannot excite a particle 
across the energy gap A so it follows that there must be 
gapless states at each end to accommodate the extra spin. 

An explicit description of the edge states requires a 
model that gives the energy bands throughout the entire 
Brillouin zone. Following Haldane [llj. we introduce a 
second neighbor tight binding model, 

ft = te L c ^ + u ^ij s if3 c L c op- ( 6 ) 

{ij)a ({ij))a/3 

The first term is the usual nearest neighbor hopping term. 
The second term connects second neighbors with a spin 
dependent amplitude, v^j — —fji — ±1, depending on 
the orientation of the two nearest neighbor bonds di and 
d2 the electron traverses in going from site j to i. = 
±1 if the electron makes a left (right) turn to get to the 
second bond. The spin dependent term can be written in 
a coordinate independent representation as z(di x d2) • s. 
At low energy (6) reduces to (2,3) with A so = 3^/3*2- 

The edge states can be seen by solving (7) in a strip ge- 
ometry. Fig. 1 shows the one dimensional energy bands 
for a strip where the edges are along the zig-zag direction 
in the graphene plane. The bulk bandgaps at the one di- 
mensional projections of the K and K' points are clearly 
seen. In addition two bands traverse the gap, connecting 
the K and K' points. These bands are localized at the 
edges of the strip, and each band has degenerate copies 
for each edge. The edge states are not chiral since each 
edge has states which propagate in both directions. How- 
ever, as illustrated in Fig. 2 the edge states are "spin fil- 
tered" in the sense that electrons with opposite spin prop- 
agate in opposite directions. Similar edge states occur for 
armchair edges, though in that case the ID projections 
of K and K' are both at k = 0. It is interesting to note 
that for zig-zag edges the edge states persist for A so — > 0, 
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FIG. 1: (a) One dimensional energy bands for a strip of 
graphene (shown in inset) modeled by (7) with t-zjt = .03. 
The bands crossing the gap are spin filtered edge states. 



where they become perfectly flatjlfj. This leads to an en- 
hanced density of states at the Fermi energy associated 
with zig-zag edges. This has been recently seen in scan- 
ning tunneling spectroscopy of graphite surfaces |l7|. 

We have also considered a nearest neighbor Rashba 
term, of the form iz ■ (s Qj a x d)c| c^g. This violates 
the conservation of s z , so that the Laughlin argument 
no longer applies. Nonetheless, we find that the gapless 
edge states remain, provided < A so , so that the bulk 
bandgap remains intact. The crossing of the edge states 
at the Brillouin zone boundary k x = it / 'a in Fig. 1 (or at 
k = for the armchair edge) is protected by time reversal 
symmetry. The two states at k x — n/a form a Kramers 
doublet whose degeneracy cannot be lifted by any time 
reversal symmetric perturbation. Moreover, the degen- 
erate states at k x = Tt/a±q are a Kramers doublet. This 
means that elastic backscattering from a random poten- 
tial is forbidden. More generally, scattering from a region 
of disorder can be characterized by a 2 x 2 unitary S- 
matrix which relates the incoming and outgoing states: 
^out = S^in, where $ is a two component spinor con- 
sisting of the left and right moving edge states 4>l], 4>ri- 
Under time reversal $i n ,out — * s y^* nt in- Time reversal 
symmetry therefore imposes the constraint S = s y S T s y , 
which rules out any off diagonal elements. 

Electron interactions can lead to backscattering. For 
instance, the term uip^dxip^ipnidxipfn, does not vio- 
late time reversal, and will be present in an interacting 
Hamiltonian. For weak interactions this term is irrele- 
vant under the renormalization group, since its scaling 
dimension is A = 4. It thus will not lead to an energy 
gap or to localization. Nonetheless, it allows inelastic 
backscattering. To leading order in u it gives a finite 
conductivity of the edge states, which diverges at low 
temperature as u _2 T 3_2A 0|- Since elastic backscatter- 
ing is prevented by time reversal there are no relevant 
backscattering processes for weak interactions. This sta- 
bility against interactions and disorder distinguishes the 
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FIG. 2: Schematic diagrams showing (a) two terminal and 
(b) four terminal measurement geometries. In (a) a charge 
current / = (2e 2 /h)V flows into the right lead. In (b) a spin 
current I s = (e/4:n)V flows into the right lead. The diagrams 
to the right indicate the population of the edge states. 



spin filtered edge states from ordinary one dimensional 
wires, which are localized by weak disorder. 

A parallel magnetic field H\\ breaks time reversal and 
leads to an avoided crossing of the edge states. H\\ also 
reduces the symmetry, allowing terms in the Hamiltonian 
which provide a continuously gapped path connecting the 
states generated by <j z t z s z and o~ z . Thus in addition to 
gapping the edge states H\\ eliminates the topological dis- 
tinction between the QSH phase and a simple insulator. 

The spin filtered edge states have important conse- 
quences for both the transport of charge and spin. In 
the limit of low temperature we may ignore the inelastic 
backscattering processes, and describe the ballistic trans- 
port in the edge states within a Landauer-Biittiker [l^ 
framework. For a two terminal geometry (Fig. 2a), 
we predict a ballistic two terminal charge conductance 
G = 2e 2 /h. For the spin filtered edge states the edge 
current density is related to the spin density, since both 
depend on nni — tie,;. Thus the charge current is ac- 
companied by spin accumulation at the edges. The in- 
terplay between charge and spin can be probed in a mul- 
titerminal device. Define the multiterminal spin conduc- 
tance by Zf = JX Gfj Vj . Time reversal symmetry re- 
quires Gft = —Glp and from Fig. 2b it is clear that 
Gf 3 = ±e/4-7r for adjacent contacts i and j. In the four 
terminal geometry of Fig. 2b a spin current I s = eV/ 47T 
flows into the right contact. This geometry can also be 
used to measure a spin current. A spin current incident 
from the left (injected, for instance with a ferromagnetic 
contact) will be split, with the up (down) spins trans- 
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FIG. 3: Feynman diagram describing the renormalization of 
the SO potential by the Coulomb interaction. The solid line 
represents the electron propagator and the wavy line is the 
Coulomb interaction. 



ported to the top (bottom contacts), generating a mea- 
surable spin-Hall voltage. 

The magnitude of A so may be estimated by treating 
the microsopic SO interaction 

h 



V so = 



T s ■ (W x p) 



(7) 



in first order degenerate perturbation theory. We thus 
evaluate the expectation value of (8) in the basis of states 
given in (1) treating ^(r) as a constant. A full evalua- 
tion depends on the detailed form of the Bloch functions. 
However a simple estimate can be made in the "first 
star" approximation: U( KK ^^ AB ^(r) — V^ p exp[iK p • 

(r — d)]/-\/3. Here K p are the crystal momenta at the 
three corners of the Brillouin zone equivalent to K or 
K' , and d is the a basis vector from a hexagon center 
to an A or B sublattice site. We find that the ma- 
trix elements have precisely the structure (3), and us- 
ing the Coulomb interaction V(r) — e 2 /r we estimate 
2A so = 4ir 2 e 2 h 2 /(3m 2 c 2 a 3 ) ~ 2.4° if. This is a crude es- 
timate, but it is comparable to the SO splittings quoted 
in the graphite literature^. 

The Rashba interaction due to a perpendicular elec- 
tric field E z may be estimated as Xr = TiVFeE z /(4mc 2 ). 
For E z ~ 50V/300nm[| this gives X R ~ .5mK. This is 
smaller than A so because E z is weaker than the atomic 
scale field. The Rashba term due to interaction with a 
substrate is more difficult to estimate, though since it is 
presumably a weak Van der Waals interaction, this too 
can be expected to be smaller than A so . 

This estimate of A so ignores the effect of electron elec- 
tron interactions. The long range Coulomb interaction 
may substantially increase the energy gap. To leading 
order the SO potential is renormalized by the diagram 
shown in Fig. 3, which physically represents the interac- 
tion of electrons with the exchange potential induced by 
A so . This is similar in spirit to the gap renormalizations 
in ID Luttinger liquids and leads to a logarithmically 
divergent correction to A so . The divergence is due to 
the long range 1/r Coulomb interaction, which persists 
in graphene even accounting for screening [2(]]. The di- 
vergent corrections to A so as well as similar corrections 
to Hvtp can be summed using the renormalization group 
(RG)[2Cj. Introducing the dimensionless Coulomb inter- 
action g = e 2 /Tivf we integrate out the high energy de- 
grees of freedom with energy between A and Ae~ e . To 



leading order in g the RG flow equations are 

dg/dl = -<? 2 /4; dA so /M = gA so /2. (8) 

These equations can be integrated, and at energy scale e, 
A so (e) = A° [l + (s°/4) log(A°/e)] 2 . Here g° and A° Q are 
the interactions at cutoff scale A . The renormalized gap 
is determined by A^ ~ A so (A^). Using an effective 
interaction g° = .74|2lj and A — 2 eV this leads to 
2A«~15°K. 

In summary, we have shown that the ground state of a 
single plane of graphene exhibits a QSH effect, and has a 
non trivial topological order that is robust against small 
perturbations. The QSH phase should be observable by 
studying low temperature charge transport and spin in- 
jection in samples of graphene with sufficient size and 
purity to allow the bulk energy gap to manifest itself. 
It would also be of interest to find other materials with 
stronger SO coupling which exhibit this effect, as well as 
possible three dimensional generalizations. 
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